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A mathematical model is developed to study the behaviour of an oscillating wave energy 
converter in a channel. During recent laboratory tests in a wave tank, peaks in the hydro- 
dynamic actions on the converter occurred at certain frequencies of the incident waves. 
This resonant mechanism is known to be generated by the transverse sloshing modes of 
the channel. Here the influence of the channel sloshing modes on the performance of the 
device is further investigated. Within the framework of a linear inviscid potential-flow 
theory, application of the Green theorem yields a hypersingular integral equation for the 
velocity potential in the fluid domain. The solution is found in terms of a fast-converging 
series of Chebyshev polynomials of the second kind. The physical behaviour of the sys- 
tem is then analysed, showing sensitivity of the resonant sloshing modes to the geometry 
of the device, that concurs in increasing the maximum efficiency. Analytical results are 
validated with available numerical and experimental data. 
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1. Introduction 

Renewed interest in wave power extraction has arisen in the recent years, following the 
need for clean energy resources alternative to fossil fuels. As a consequence, a lot of work 
has been done recentl y to d i scove r effective design layouts for wave energy converters 
(see iMcCormickl 120071 : ICrud 120081, for an extensive review of wave energy conversion 
techniques). The work of lFollev et~al\ (|2005l ) and IWhittaker et all (|2007h on a number 
of different devices showed that an effective solution to harness energy from water waves 
is represented by oscillating wave energy converters. Their simplest form is that of a 
buoyant flap hinged at the bottom of the ocean and has been named Oscillating Wave 
Surge Converter (OWSC). Under the action of incoming incident waves, the flap is forced 
to a pitching motion about the hinge, that is eventually converted into energy with 
the use of a generator. The converter is naturally placed in shallow waters, where the 
pitching of the flap combines well with the amplified horizontal motion (surge) of the flui d 
particles (hen ce t he name OWSC) . During the wave tank tests of lWhittaker et all (|20071 ). 
Henry (2008) and Ivan 't Hofj (l2009h . peaks m the hydrodynamic actions on the converter 
occurred at certain frequencies of the incident waves. This resonant phenomenon is known 
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to be due to the excitation of the transverse sloshing modes of the channel (jShemer et 



19871 ). However, there is still uncertain ty on how the re sonant transverse waves interact 



with the device and affect its efficiency ( van't H off 2009). Here we address this unresolved 
issue by developing a three-dimensional mathematical model of an oscillating wave energy 
converter, namely the OWSC, inside a channel. Although much analyti cal modelling 
has been done f or a large number of wave energy converters (WECs) (see IFalnesI 2002 ; 
Mei et al. 20051 iMcCormick 2007 . for a collection of theories on WECs), few analytical 



studies are available in the literature concerning OWSCs. 

Indeed, the hydrodyn amics of flaps has been w idely studied in the past due to their 
usage as w avemaker s fseelDean fc Dalrvmplell99ll for a compendium of wavemaker theo- 
ries), with lMilgraml (| 19701 ) being among the first to suggest the employment of a moving 
paddle not only as a wave generator, but also as a wave energy a bsorber. Later on, 
the hydrodynamic characteristics of flat pl ates have been s t udied by IParsons fc Martin 
( 19921 1 1994 , Il995l ) and more recently by lEvans fc Porter ] dl996l) for both subm erged 
and surface-piercing plates. In particular, the work of IParsons &: Martini (|l992l ) is a 
seminal example of the application of hypersingular integral equations to problems in 
water waves. Nevertheless, all these analytical studies deal with two-dimensional (2D) 
geometries. Hence they are inadequate to describe the complex three-dimensional (3D) 
behaviour of the wave field around the flap noticed in the experi mental ana l ysis o f the 
OWSC. A non- linear 3 D analytical model has been developed bv lMei et ali (|1994f) and 
Sammarco et all (Il997l) fo r the flap gates of the barrier s designed to protect Venice from 
flooding. Mei et g/] ~( 1994h and lSammarco et al. ( 1997 ) showed that trapped modes ex- 
ist near a barrier of oscillating flaps in a channel. Since the narrowly spaced flaps span 
the whole channel width, this s ystem is closed and can resonate only subharmonically 
(i.e. non- linearly). The work of ISammarco et all (|1997| ) aims at reducing the resonant 



effects that could lead to undesired large oscillations of the flaps, thus undermining the 
effectiveness of the barrier. 

Going back to wave energy, a wise design of the OWSC should exploit the bene- 
ficial effect of resonance in enhancing the amplitude of oscillation of the flap, thus 
allo wing to capture more wave p ower. For bodies of va r ious geometries i n a c han- 
nel. lEatock Taylor fc Hrmel (|l985l) . lYeung fc Sphaierl (|l989l ). lLinton fc Evansl (|l992l ) and 
IChenl (| 1994T) already showed that linear resonance mechanisms influence the radiation 
and diffraction properties significantly. This influence manifests itself by the presence of 
peaks in the hydrodynamic characteristics of the system, at periods corresponding to the 
resonance of symmetric transverse standing waves within the channel lateral walls. As 
already mentioned, similar dynamics has been shown to occur also during wave tank tests 
of the OWSC, rising the following questions: (i) how linear resonance of transverse waves 
is triggered for a typical OWSC configuration in a channel, (ii) how do the transverse 
sloshing modes interact with the OWSC and modify its performance, and (iii) to what 
extent the results obtained in an experimental wave tank can be applied to predict the 
behaviour of the device in the open ocean. 

A mathematical model of the OWSC in a channel is developed in the following sections 
to address the above points. The fluid is assumed to be inviscid and incompressible, the 
flow irrotational, the perturbation time harmonic. In <|2] the governing equations are 
scaled and linearised according to first-order potential flow theory. The hydrodynamic 
problem is then decomposed into a radiation and a diffraction problem, which are solved 
separately using the hypothesis of thin plate and eventually summed up to obtain the 
total potential. The solution is achieved with the integration of a hypersingular problem, 
resulting from the application of the Green theorem to the physical domain. In Sj3] the 
hydrodynamic characteristics of the system, such as the excitation torque, added inertia 
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Figure 1. Geometry of the converter and the channel: (a) section, (6) plan view. 

torque and radiation damping are fully investigated and the efficiency of the OWSC is 
assessed. Analytical results agree well with available numerical and experimental data. 
It is shown that all the above quantities are characterised by the occurrence of spikes 
at the cut-off frequencies of the transverse sloshing modes. The latter are excited in the 
channel by the diffraction of the incident wave field at the edges of the flap and influence 
the performance of the OWSC. Finally, a parametric analysis reveals how the transverse 
modes can increase the efficiency of the device, depending on the geometry of the system. 



2. Analytical model 

2.1. Governing equations 

Referring to figure QJa, b), consider an OWSC in water of depth h', where the prime 
indicates a physical dimensional quantity. For the sake of generality, we shall represent 
the converter by a rectangular flap of width w' and thickness 2a', hinged along a straight 
axis upon a rigid platform, at a distance d from the bottom of the ocean (see figured]!, 6). 
The device is located in the middle of a straight channel, whose impermeable walls are 
placed at a mutual distance b' and extend to infinity to either side (see figureQJ). Incident 
waves are coming from the right with wave crests parallel to the OWSC and set the flap 
into an oscillating motion. This is eventually converted into useful energy by means of a 
generator linked to the device. Note that, due to the mirroring effect of the side walls, 
the layout of figure can also be employed to model the behaviour of an infinite 
array of identical OWSCs in the open ocean, subject to beam seas. The converters in 
this array have spatial period b' and oscillate at unison. Returning to the channel layout, 
let us now define a plane reference system of coordinates x' = (x',y',z'). Let x' lie on 
the centre line of the channel, y' along the axis of the OWSC at rest position and z' 
rise up from the undisturbed water level z' — 0, positive upwards. The flap is able to 
oscillate on the vertical plane (x' , z') about the horizontal axis at (x 1 , z 1 ) = (0, —hi + d), 
thus representing a system with one degree of freedom, i.e. pitch. Its time-dependent 
amplitude of rotation 9' = 6'(t') is defined positive if counter-clockwise; t' denotes time. 
We also assum e that the bottom foundation does not oppose the flap movement (ideal 
foundation, see Sammarcolll996l ) and that the hinge is frictionless. As in many gravity- 



4 



E. Renzi and F. Dias 



wave problems (see for example iMei et al\ 120051 ) the fluid is deemed to be inviscid and 



incompressible and the flow irrotational. Hence there exists a potential &'(x',y',z\t' 
for the velocity field v' = V'$', where V'(-) = {(-).x', ('),y', (")>«'} ^ s the narj la operator. 
Subscripts with commas denote differentiation with respect to the relevant variable. The 
complete set of nonlinear equations governing the system is detailed in Appendix |A"1 



Since the resonance of the channel sloshing modes is a linear phenomenon (jShemer et 



19871) , here we shall consider a linearised theory. The first-order analysis presented in this 
paper constitutes the basis for further investigation of higher-order phenomena occurring 
in the system (see $3]) . Now let A' denote the amplitude scale of the wave field and g the 
acceleration due to gravity. Then introduce the following non-dimensional variables 

(x, y, z) = (x\ y', z')/b', t = Jin, C = ('/A', $ - (y/gVA) ~* = 6'/e, (2.1) 

where e = A'/b' is the physical scale of the angular rotation, and constants 

{h,a,w) = {h',a',w')/b'. (2.2) 

The linear governing equations are derived in the limit of small-amplitude oscillations of 
the flap, i.e. e -C 1, by considering the terms 0(1) in the fully nonlinear system (|A 2[) - 
(|A 10p scaled according to (|2.1I) . The potential is governed by the Laplace equation 

V 2 $ = 0, (x,y,z)en (2.3) 

where Q is the fluid domain. On the free-surface, the kinematic-dynamic boundary con- 
dition reads 

+ *,« = <), z = 0. (2.4) 
The no-flux conditions on the solid boundaries require 

$ )3/ = 0, y = ±1/2 (2.5) 

on the channel lateral walls and 

$ 2 = 0, z = -h (2.6) 

on the bottom. Allowing only tangential motion along the lateral surfaces of the OWSC 
yields 

^ x --e t {t){z + h-c)H{z + h-c), x = ±a, -w/2 < y < w/2, (2.7) 

where H denotes the Heaviside step function. Finally the linearised equation of motion 
of the flap reads 

/0 pw/2 
/ [$ t (a, y, z, t) - $ t (-a, y, z, t)] (z + h-c) dydz + T e (t), 
-h+c J -w/2 

(2.8) 

where C = wa(h — c) 2 — S is the elastic moment of inertia due to net buoyancy, S — 
S' I (pb' 4 ) and / = I' /(pb 15 ) are respectively the first and second moment of inertia of the 
flap and T e = T^/ipgb' 4 ) is the first-order torque exerted by the generator on the flap, 
all in non-dimensional variables; p is water density. Without loss of generality, we shall 
assume % to depend on time via 9(t) and to be partly inertial, partly elastic and partly 
damping (see Mei et aZll2005l for a similar example) 



Te(t) = -fi p to <tt {t) - C p to 6{i) - v vto >t (t). (2.9) 

In the latter expression p pto and C p t are respectively the inertial and elastic characteris- 
tic of the generator, while v pto is the energy extraction rate. Now note that the physical 
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= O(A') (|Henrvll2008h . so that a = a'/V = 0(e) < 1. 



Hence the thickness of the OWSC can be ne glected if compared to t he reference length 
scale b' and the hypothesis of thin obstacle ( Linton fe Mclverll200ll) can be applied to 



evaluate the potential $. As a consequence, a = in the boundary condition on the 
lateral surfaces (12.71) and in the right-hand side of the equation of motion l|2. 



2.2. Solution of the radiation and scattering problems 

In this section solution is found to the Laplace equation (|2.3p with linearised boundary 
conditions (|2.4I) on the free-surface, (|2.5p on the channel walls, (|2.6[) on the bottom and 
(|2.7p on the flap lateral surfaces. Once having solved this system, expression (|2.8p will 
be employed to determine the motion of the flap. Assume that the flap undergoes simple 
harmonic oscillations of frequency u> — us 1 \JV~j g about its rest position, so that 

0(t) =Rc{Qe- iult } , (2.10) 

where Re is the real part and 8 the complex amplitude of rotation. Position (|2.10l) allows 
to separate the time factor in the time-dependent variables, so that the velocity potential 
can be rewritten as 

<!>(x,y,z,t) = 4>(x,y,z)e- i "\ (2.11) 

where the symbol Re is omitted for the sake of brevity. Due to the linearity of the 
problem, the sp atial potential <j) is analysed by resortin g to the classical decomposition 
(see for exa mple iLinton fc Mclverll200ll: iMei et all [20051) 

q) = (f> R + (j) s . (2.12) 

In the latter, cf> R is the solution of the radiation problem, in which the flap is set to oscil- 
late without incoming incident waves, and cj) is the solution of the scattering problem, 
where the flap is held fixed in incoming waves. In turn, the scattering potential can be 
decomposed into 

s = 7 + <^, (2.13) 

where 

tf&y^z) = '-^-^ CO sh k(z + h) C - ikx (2.14) 

us cosh kh 

is the potential of the i ncident wave an d (j> D the potential of the diffracted waves (see 



Linton fc Mclverll200ll ; IMei et aZ.I I2005T ). In (|2. 14[) A and k are respectively the non- 
dimensional incident wave amplitude and wavenumber, the latter depending on the wave 
frequency according to the dispersion relation ui 2 — ktanhkh. In the following, the 
radiation and scattering problems will be investigated and their results summed up to 
obtain the total wave field. 

The boundary-value problems for the radiation potential <fs R and the diffraction po- 
tential <j> D can be determined by employing the factorizations (|2.10[) - (|2.11[) and the 
decompositions (|2.12|) - (|2.13[) into the governing equations (|2.3j) - (|2.7j) . As a consequence, 
4> R and 4> D must satisfy the Laplace equation 

V 2 ^' D )=0, (x,y,z)en, (2.15) 

where the shorthand notation cf>( R ' D ) denotes either potential, the kinematic-dynamic 
boundary condition on the free-surface 

the no-flux condition on the bottom 



4 R > D) = o, z = -h 



(2.17) 
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and on the channel walls 

0^=0, y = ±l/2 (2.18) 
and the kinematic conditions on the lateral surfaces of the flap 

iujQ (z + h — c)H(z + h — c) 1 . n 

v v ; | , x = ±0, -w/2 < y < w/2. (2.19) 

Finally, (f>( R ' DS > must be outgoing at large |x|. Separation of variables 

oo 

^ a > D Hx,y ) z) = y £ iV ^ D \x,y)Z n (z), n = 0,l,2,..., (2.20) 







f ] 




I 42 J 


H 







yields the well-known vertical eigenmodes (see for example iMei et a/.l l2005) 

ry ,S \/2COSh K n (z + h) 

Z n {z) = ■ — y/2> "- = 0,1,2,..., (2.21) 

(h + uj~ 2 sinh K n h) 

which satisfy the orthogonality relation 

o 

Z n (z)Z m (z) dz = 5 nm , (2.22) 

l-h 

Sum being the Kronecker delta, n, m G N. In Z n (|2.2ip 

K = k, K n = ik n , 7i = l,2,... (2.23) 

are the solutions of the dispersion relationships 

uj 2 = fctanhfcfr, uj 2 — —k n tan k n h, n = l,2, ... (2.24) 

respectively. Hence an equivalent form of (|2.2ip for n > is Z n — \[2 cos fc n (^ + h)/(h — 
lj~ 2 sin 2 knh) 1 / 2 , i.e. an oscillating function of z. By making use of (|2.20p and (I2.21[) , the 
Laplace equation (|2.15[) becomes 

(V 2 + kI) ^ = 0, (2.25) 

governing the plane radiation and diffraction potentials (p n R,D ^ (x, y) in the 2D fluid do- 
main E (see figure [T]), where n — 0,1,2,... will be omitted from now on for brevity. 
Note that (12.251) is an ordinary Hclmholtz equation when n = 0, for which Kq — k 2 , and 
a modified Hclmholtz equation when n > 0, for which k 2 = — k 2 . With the same steps 
as above, the no-flux condition on the channel walls (I2.18[) gives 

<P%Z D) =0, 2/ = ±1/2- (2.26) 

Now consider the kinematic conditions on the flap (|2.19[) . Using the factorization (|2.20p . 
multiplying both sides by Z m , integrating along the vertical line from z = — h to z = 
and using the orthogonality relation (|2.22p yields 



<« If Vf n 



D \ 1 A d n 



n.x 



x = ±0, -w/2 < y < w/2. (2.27) 



In the latter expression 

V = iluO (2.28) 
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is the complex angular velocity of the flap, while 



v2 [n n (h — c) sinh n n h + cosh n n c — cosh.K n h] 



(h + uj 2 sinh 2 K n Ji) 



1/2 



and 



d = 



k(h + ui 2 sinh 2 kh) 
\pl uj cosh kh 



1/2 



d n = 0, n = 1,2, . . 



are real constants depending on the modal order n. Finally ipi^' D ^ must be outgoing dis- 
turbances along the channel. In summary, the 2D n-th mode potentials tp^'^ must solve 
the Helmholtz equation (|2.25[) with the boundary conditions (I2.26p - (|2.27[) and be outgo- 
ing. Here we shall solve this boundary- value problem by employing an integral-equation 
method based on the appli cation of Green's t heore m wit h an appropriate chann el Green 
function. In the literature lAchenbach fe Lil (Il986h and IMartin fe Ri^ol (Il989b already 
used a simi l ar pr ocedure to solve plane sound-wave problems, while Parsons fc Martini 
( 19921 1199411 19951) applied this method to the scattering and trapping of surface waves by 
plates. Later, IMartin fc Farinal (|1997j) used the integral-equation approach to solve the 
radiation problem for a heaving submerged horizontal disc, but without bounding walls 
and in water of infinite depth. In this paper, the boundary-value problem (|2.25l) - (|2.27p 
is solved in Appendix [B] by applying an integral-equation technique based on the decom- 
position of the Green fu nction into a singular part and an analytical remainder. Such a 
procedure is inspired bv lChen ( 1994 )'s analysis of the tank Green function for numerical 
models. The integral-equation technique allows to solve the plane boundary-value prob- 
lems for tpi R ' D ' > and ultimately to express the potentials ^ R ' D ^ (|2.20p in an appealing 
fast-converging semi-analytical form (see Appendix [B] for details). The radiation potential 
is given by 



(2.29) 



(2.30) 



iwV 



M 



-t-oo 



E / i 1 



n=0 



p=0 

H 
x ■ 



m— — oo 



,2\l/2 



U 2p (u) 



(^K n y^x 2 + (y — iwu — m) 2 ^ 
\/ x 2 + (y — ±wu — m) 2 



du, (2.31) 



where is the Hankel function of the first kind and first order and the U2 P are the 
Chebyshev polynomials of the second kind and order 2p,p = 0, 1 . . .M G N. Finally, the 
0i (2p)n are the complex solutions of the linear system of equations (|B 231) which ensures 
that the boundary condition on the plate (|2.27|) is satisfied. This system is solved numer- 
ically with a collocation scheme (see again Appendix [Bj , therefore the solution is partly 
numerical. Note that the radiation potential (|2.31j) is an odd function of x, as suggested 
by the antisymmetric pitching motion of the flap. Physically, (12.311) represents the sum of 
outgoing waves within the waveguides imposed by the lateral walls and is similar in form 
to the soluti ons of other problems of water waves p ropagating past sharp obstacles (see 
for example Evans fc Porterlll997t lLovas et aZ.I l2010l who also employ Chebyshev poly- 
nomials) . The inner series in (|2.3ip represents a sum of disturbances radiated by sources 
arranged periodically along the y axis and converges as m -3 / 2 for large m. Numerically, 
summation from m — —20 up to m = 20 ensures convergence of the series with relative 
error O(10~ 3 ) close to resonance and O(10 -4 ) elsewhere. The central series in (|2.31l) is 
truncated at M < oo and its coefficients are calculated numerically as described in Ap- 
pendix |5] Numerical investigation showed that M ^ 4 is sufficient to ensure convergence 
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with virtually no relative error. The outer series in (|2.31|) over the vertical eigenmodcs 
was found to converge quickly already for n = 4 with a relative error O(10 -5 ). Hence the 
method of solution adopted here, based on the appropriate decomposition of the Green 
function (|B 4[) - (|B 6[) . proves to be reliable and computationally efficient. Turning to the 
scattering problem, the diffraction potential is given by 

M +00 



p—0 m— — OG 1 

{k^/x 2 + {y- iron - m) 2 ) 



x 2 + (y — iwti — m) 2 



du. (2.32) 



In the latter expression the /3(2 P )o are the complex solutions of the linear system of 
equations (|B 23|) for n = 0, which ensures that the no-flux condition on the flap (|2 . 19[) 
is satisfied. Thanks to linearity, the total potential = cf> R + cf) D + (j) 1 is obtained by 
simply summing up the solutions (f> R (|2.31|) and D (|2.32j) of the radiation and diffraction 
problem, respectively, to the incident wave potential (j) 1 (|2.f 4[) . The equation of motion 
of the flap can now be solved. 

2.3. Flap motion 

2.3.1. Hydrodynamic parameters 

Consider the law of motion of the flap (|2.8|) . By using the factorizations (|2.10|) and 
(|2.11[) for 8 and $ respectively, applying the decomposition cf> = cf> R + <p D + cf) 1 , employing 
the solutions (|2.31[) and (|2.32[) for <j) R and cf) D respectively and expressing % with (|2.9p . 
d2~8]) yields 

[-w 2 (7 + /i + ^ pto ) + C + C p to - iu(v + u pto )] Q = F. (2.33) 

In the latter expression Q is the unknown complex amplitude of rotation of the flap, 
while 




/i = — He ; \ ; (2.3J) 

is the added torque due to inertia and 



n=Q 



UJ1TW 



+00 



Im <^ Y, a 0nfn \ (2.35) 



^n=0 



is the radiation damping, both depending on the solutions ao n to the linear system (|B 23|) 
for the radiation problem. Finally 

F = -—iojA (3 00 f (2.36) 

is the complex excitation torque, corresponding to the action on the plate as if it was 
held fixed in incoming waves. 

2.3.2. Extracted power and capture factor 

Expression (|2 .33[) formally describes a simple harmonic oscillator in the frequency 
domain; hence the average extracted power from the OWSC over a period T = 2tt/uj is 

1 r\„ OA*.*- 1 ^ Pto \F\ 2 



p =r (v P toO,t) 6, t dt = ' 5 — ' ( 2 - 37 ) 

1 Jo 2 [C + Cpto - (I + /J. + Mpto) w 2 ] + (y + i/ pto ) u 2 
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where again 9 = RejOe -1 " 4 } and is the solution of (|2.33|) . Returning to physical 
variables via (|2.1|) . the power take-off expressed in Watts (W) is then 

P' = P A' 2 b' 3/2 g 3/2 P [W] . (2.38) 

The rate of power take-off (|2.37p can be maximized according to the optimizing criteria 



C + Cpto 



l) = Uq = a — ■ , v pto = v. (2.39a, b) 

y I + /i + \i vto 

In the latter, condition (12.39b ) corresponds to the tuning between the incident wave 
frequency oj and the natural pitching frequency coq of the flap, while (|2.39b ) requires the 
equality of the rate of power take-off v p to to the radiation damping v (|2.35p . A quick 
numerical assessment of the tuning condition (|2.39b ) for a typical OWSC configuration 
(see §3.ip reveals that the natural pitching period Tq = 2-tt /ujgy/b ' / 'q ~ 20 s is larger 
than the ordinary wave periods on site, about 5 — 15 s (jHenrv 2008). Hence the tuning 



condition (|2.39b ) is not achieved in normal operating circumstances. Away from body 
resonance, the extracted power (|2.37|) can still be optimized by adjusting the damping 
rate v pto such as dPjdv pto — 0, which yields 



^pto — ^opt 



[C + C p to -(I + fi + Upto)^ 



(2.40) 



as the optimum damping rate (sec Falncs 2002). Then the optimum power P opt is obtained 



by substituting v pto — v opt in (12.37[) . Now, the rate of extracted power and the optimum 
power are functions of the wave climate, so that the amount of energy extracted from a 
given floating body strongly depends on the sea state. Therefore both the rate of power 
take-off and the optimum power are not reliable indicators of the efficiency of the system 
(|Crud 2008). A more rational efficiency assessment is done by considering the capture 



factor, defined as the ratio between the optimum power extracted per unit flap width 
and the power of the incident wave field per unit crest length, i.e. 



P' 

op 



where A' a = A A' and 



C' = — (l + ^ ) [m/s] (2.42) 
9 2k' \ smh2k'h' ' 



are respectively the amplitude and the group velocity of the incident waves in physical 
variables. In the following, all the parameters describing the performance of the OWSC 
will be assessed in physical variables. 

In the next sections, the analytical theory will be validated against numerical and 
experimental models. Discussion will follow on the physical behaviour of the system and 
further analysis will investigate the dependence of the system on its main parameters, 
such as the frequency of the incident waves and the width of the flap. Resonance will be 
shown to occur at special frequencies, namely the cut-off frequencies of sloshing waves 
in open channels. The response curves of the system dynamic characteristics show that 
the OWSC performance is largest under resonance conditions. 
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Figure 2. Behaviour of (a) added inertia torque and (6) excitation torque with wave period in 
physical variables. The solid line refers t o the analytical values of i]2,3.1l while the markers refer 
to the numerical solution of Ivan 't Hofl (pOOOh for various thicknesses of the flap. 



3. Discussion 

3.1. Numerical and experimental comparison 

In ord er to validate the theory, comparison is made with the numerical results of lvan't Hofl 
(j2009h . In the analytical model, the width of the flap is w' = 18 m, while the width of 
the channel is b' = 91.6 m, giving a blockage ratio w = w'/b' ~ 0.2. The depth of the 
channel is h! — 10.9 m and the distance between the bottom and the hinge is c' = 1.5 m. 
The amplitude of the incident wave is A' = 0.3 m, which gives e = A' /b' ~ 0.003 <C 1. 
This justifies the assumption of linear theory for this case. The nu merical data have been 
obtained with the software package WAMIT (|WAMIT Inc.ll2008h for the linear analysis 
of the interaction of surface waves with offshore structures. The geometry of the flap 
used in the numerical simulations is that o f an 18 m wide r ectangular box standing on 
a triangular prism pointing downwards fsee lvan't Hofl 2009, for specifications). The flap 
in turn is hinged on a 1.5 m high rectangular platform lying on the bottom of a 10.9 m 
deep ocean. In the numerical calculations of Ivan 't Hod (2009) the device is placed in 
open water, so that the model does not simulate the effect of the lateral walls. Figures 
Ufa) and[U&), respectively, show the values of the added inertia torque fjt' = pb' 5 fi and 
the absolute value of the excitation torque \F'\ — pgA'b' 3 \F\ in physical variables versus 
the wave period T", for both the analytical and the numerical model. While in the an- 
alytical model the flap has zero thickness, in the numerical model the simulations have 
been carried on with different non-zero values of the total thickness 2a' (see again figure 
[2]z,&). The agreement between both models is satisfactory. However for the added inertia 
torque p! (figure [2ji) , the analytical values are slightly larger than the numerical ones. 
This is likely due to the differences in the shape of the flap, a perfect rectangular box in 
the analytical model and a rectangular box sitting on a triangular base in the numerical 
model. Also, note that the analytical values are larger than the numerical ones next to 
the peaks of either the added inertia torque (T 1 ~ 9.6 s in figure HJi) and the excitation 
torque (T" ~ 5.7 s in figure W?)- This is a co nsequence of the bounding effect of the 
lateral walls, absent in the numerical model of lvan't Hoi d 2009), that will be discussed 
in depth in £13.21 Finally, the numerical results in figure [2ja,6) clearly show that the 
hydrodynamic coefficients do not vary noticeably while varying the flap thickness. The 
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Figure 3. Behaviour of the excitation torque with the period of the incident waves in physical 
variables. The squares refer to the experimental results of Henrv (2008) for the 20 th scale model of 
the Oyster™ WEC, while the solid line refers to the equivalent analytical model as developed in 
section Sj2] Discrepancy between data at T ~ 10 s can be explained by the end-to-end resonance 
of long-crested waves in the wave tank. 




initial assumption of neglecting the thickness of the flap to calculate the hydrodynamic 
actions, namely the thin-plate approximation, is then valida ted. Th e theo retical results 
of <|2] are further validated by comparison with the data of iHenrvl (2008). These have 
been obtained during an experimental campaign at Queen's University Belfast on a 20 th 
scale model of the Oyster™ wave energy converter developed by Aquamarine Power Ltd 



Q (see lFollev e7aDl2007t Iwhittaker et olll2007t iHenrv et aZjEnoh having similar dimen- 



sions to those already specified for the theoretical models. In figure[3]the behaviour of the 
excitation torque \F'\ with the incident wave period T" is reported for both the equiva- 
lent analytical box-shaped flap and the experimental model of the Oyster™ WEC. Very 
good agreement is found between the two sets of data: the theoretical model described 
in $2] is also capable to simulate the behaviour of a real - and more complex - OWSC. 



3.2. Resonance 

After having shown a good agreement between the analytical and the numerical/exper- 
imental data, we shall investigate in more depth the 3D dynamics of the flap in the 
channel and the influence of the channel lateral walls in enhancing the efficiency of the 
converter. First, let us consider the diffraction problem and analyse the modifications 
induced by the 3D wave field in the 3D channel model of $2] with respect to a simplified 
reference model, in which b' = w , i.e. the flap spans the entire width of the channel. Here 
the behaviour of the fluid is clearly two dimensional, so that thi s layout will be referred 



to as the 2D channel. As already demonstrated in the literature (jLinton fc Mclverll2001 



Fames! 2002 ). the excitation torque acting on the flap in the 2D channel is, in physical 
variables, 



F 2d 



2pgA' w' 



h! - c' 



tanh k'h' 



cosh fc'i 



cosh k'h' 



k' 2 coshk'h' 



[Nml 



(3.1) 



f www.aquamarinepower.com 
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Figure 4. Behaviour of the excitation torque with the period of the incident waves in physical 
variables. The solid line refers to the 3D value F', calculated for the geometry of £13.11 . The 
dashed line refers to the 2D value F^d (|3.1[1 for the simplified 2D channel, where b' — w' = 18 m. 
The dash-dotted line represents the 2D long- wave limit (|3.2p . Note that the 3D torque has a 
maximum at T' = 5.7 s and spikes at T' = 4.5 s and T' = 9.6 s. 



In figure@]the excitation torque is plotted against the wave period T 1 in physical variables, 
for either the 3D-channel model of Sj2] and the 2D-channel reference model. While the 
2D excitation torque regularly increases with the period up to the asymptotic long-wave 
value 



pgw'A' (ti - c'Y [Nm] ~ 4.7 ■ 10 6 Nm, 



(3.2) 



the 3D torque \F'\ reaches its maximum value max {|F'|} ~ 4. 75 TO 6 Nm at a finite period 
T' = 5.7 s. Then it decreases while moving towards longer periods, with a secondary spike 
at T' = 9.6 s. Hence an immediately noticeable effect of the 3D dynamics is the lowering 
of the peak period to finite values. Though the two curves of figureS]have a very different 
aspect, they show the interesting property that max {\F\ } / max {F^} — 0(1) (see again 
figure 0]). Hence to further study the 3D effects on the system, it is convenient to consider 
the non-dimensional ratio 



(3.3) 



between the excitation torque in the 3D channel and the maximum excitation torque 
in the 2D channel. In figure [5ja) / is plotted versus the non-dimensional wavelength 
A = A' /b' of the incident wave, for the same geometry as ab ove. The graph o f figure 
[5ja) exhibits the typical spiky behaviour already noticed by lAchenbach fc Lil (1986) 
when studying the reflection coefficients for normal inci dence of sound waves on an ar - 
ray of screens. Similar s pikes have also been fo und bv lEatock Tavlor fc Hungl (J1985J), 
lYeung fc Sphaierl (|l989l ). iLinton fc Evansl (|l992t ) and[Chen| dl99^) in the analysis of the 
hydrodynamic coefficients of bodies of various shapes in a channel. Note that in figure 
[5fa) the spikes, either in the form of local maxima or changes in slope, occur at special 
values of A, i.e. A = 1, 1/2, 1/3, . . . The latter correspond to the well-known values of 
the cut-off wavelengths of the symmetric sloshing modes in a channel (see for example 





Figure 5. (a) Non-dimensional excitation torque ratio / (|3.3[) versus the wavelength of the 
incident wave for the geometry of £|3 . 1 1 - Note the spiky shape of the curves, with spikes occurring 
at the resonant frequencies Ai = 1, A2 = 0.5 and A3 = 0.333 . . . marked by vertical dash-dotted 
lines, (b) Group velocity ratio c gm (|3.5[l versus wavelength for the first three sloshing modes. Note 
that the group velocity of each sloshing mode decays as soon as the relevant cut-off wavelength 
A m (|3.4[) is approached. 
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The sloshing modes are transverse standing waves whose crests are parallel to the axis 
of the channel, oscillating back and forth between the two bounding lateral walls. The 
physical behaviour of the sloshing waves is strongly depende nt on the w i dth of the channel 
and on the wavelength of the incident wave. As shown bv IMei et al] (|2005l ). if X'/b' = 
A < A m , then the m-th sloshing component propagates along the channel, carrying the 
associated wave energy towards infinity. On the other hand, if A > A m , the amplitude 
of the m-th sloshing mode decays exponentially while moving along the channel. In 
our case, this means that the energy associated to these modes is trapped near the 
OWSC and available for extraction. A further insight on this fundamental influence of 
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the sloshing modes can be achieved by consider ing the ratio bet ween the TO-th sloshing 
mode group velocity C' = C'yl — (mA) 2 fsee lMei et a/.ll2005l ) and the group velocity 
of the incoming waves C (|2.42l) . i.e. 



c gm = C' g JC' g = y/l - (mA) 2 . (3.5) 

Recall that the group velocity is the speed of energy transport. As soon as A — > A m = 1 /m, 
c gm and the energy of the m-th sloshing mode is no longer propagating along the 
channel: the sloshing mode is trapped near the flap. Figure EJo) shows the behaviour of 
the group velocity ratio c gm (13.5[) for the first three sloshing modes versus the wavelength 
of the incident wave. Let us first consider the group velocity ratio for the mode m = 3 
(dotted line). When A — > A3 = 1/3, the group velocity ratio c 9 3 quickly decays to zero, 
so that the third sloshing mode is now trapped near the flap. The trapping of the energy 
associated to this mode determines a sudden increase of the excitation torque acting 
on the flap, which results in / having a spike at A = 1/3, as shown in figure Eta). By 
further increasing A, the group velocity ratio of the second sloshing mode quickly decays 
too, being c g2 = for A = A 2 = 1/2 (dashed line of figure [5f>). Due to the trapping 
of the second sloshing mode, the excitation torque acting on the flap further increases, 
leading to the peak of / occurring exactly at the cut-off frequency of the second sloshing 
mode, A = A2 (see figure |5Ji). Additionally increasing A makes the incoming waves so 
long that they pass through the flap almost as a uniform swelling. As a consequence, the 
net pressure acting on the flap diminishes, leading to a decrease of the excitation torque 
F which also reflects on / (see again figure Eb). Nevertheless, as soon as A Ai = 1, 
c gi ~> (see the solid line in figure [S^) and the first sloshing mode is also trapped near 
the flap. This results in / having a ultimate spike at A = 1, as shown in figure [SJa). 
The physical picture is now clear. The incoming plane waves, initially 2D, impact the 
flap in the middle of the channel. Differently from the 2D channel, in which the incident 
waves are totally reflected back to the source, in the 3D channel the waves are also 
transmitted beyond the flap, due to the lateral gaps between the body and the channel 
walls. Diffraction at the edges of the flap modifies the behaviour of the wave field and 
dominates the shape of the hydrodynamic coefficient curves, making them different from 
the 2D scenario (see figure 2}. As a consequence of this enriched dynamics, transverse 
sloshing modes, a feature of the channel, are ultimately excited. Hence the initial 2D 
motion shatters into a series of symmetric 3D sloshing waves, resonating each at its own 
cut-off frequency. When a sloshing mode resonates, the relevant energy is trapped near 
the flap, resulting in an increase of the excitation torque. Similar conclusions can be 
drawn also for the radiation problem, where the waves are generated by the pitching 
motion of the plate. Figure H^a, 6) shows the plots of the added inertia torque fi (|2.34|) 
and the radiation damping v (|2.35p versus the wavelength of the generated waves in 
non-dimensional variables. Note that spikes occur here at the same resonant wavelengths 
as those of / (|3.3p for the scattering problem (see again figure [SJ . This correspondence 
is due to the relationship between the radiation and scattering problems 



mw . Vk (2kh + smh2kh) 



1/2 



8 coshfc/i 00 

shown in AppendixfCj for which the excitation torque acting on the fixed flap in incoming 
waves is associated to the radiation potential of the pitching motion of the body. We now 
investigate the OWSC performance. 




Figure 6. Behaviour of (a) added inertia torque and (b) radiation damping with the wave- 
length of the generated waves in non-dimensional variables for the geometry of £|3 . 1 1 Vertical 
dash-dotted lines indicate the resonant frequencies of the transverse sloshing modes. 
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Figure 7. Behaviour of (a) amplitude of rotation and (b) capture factor with the period of 
the incident wave in physical variables for the geometry of £|3. II In panel (b) the horizontal 
dash-dotted line indicates the average capture factor. 



3.3. Wave power extraction 

After having analysed the radiation and scattering problems separately, we now turn to 
study the performance of the OWSC in terms of wave power extraction. First, let us 
consider the flap pitching amplitude 9, solution of (|2.33p . For easiness of representation, 
but without loss of generality, in expression (|2.33l) we take the generator inertial and 
linear terms /K pto = C pto = and make use of the optimizing criterion (|2.40l) for the 
rate of power take-off, such that v vto = v Q pt- Finally, we use the values of the moment of 
inertia I and buoyancy factor C for a typical configuration of the Oyster™ wave power 
device (private communication with Aquamarine Power Ltd). In figure [7ja) the absolute 
value of the flap pitching amplitude is plotted versus the period of the incident wave in 
physical variables. Note that the amplitude of oscillation of the flap increases smoothly 
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while moving to longer waves, not showing the characteristic spiky behaviour of the 
hydrodynamic coefficients. Since O is related to the displacement in the x direction via 
the pitching movement of the flap, it seems unaffected by the sloshing modes resonating 
in the orthogonal y direction. Let us now consider the capture factor Cp, defined by 
(|2.4ip as a measure of the OWSC efficiency in catching the energy available in incident 
waves. In figure the capture factor is plotted against the incident wave period. 
Comparison between figures 2] and reveals that both the excitation torque and the 
capture factor have a peak at T' = 5.7 s, which corresponds to the resonant period of the 
second sloshing mode. Hence the capture factor is dominated by the excitation torque, 
which is the power generating action. Overall, figure WU>) reveals that the average capture 
factor of the 3D-channel model is larger than the maxim um theoretical value Cf = 5 of 
the 2D-channel model (see for example Mei et alll2005l ). This shows again that the 3D 



dynamics described in §3.21 altogether increases the efficiency of the device with respect 
to the 2D layout. Finally, the maximum capture factor at resonance is maxCj? ~ 0.80, 
indicating that the interference between the transverse sloshing modes occurring in the 3D 
channel can further enhance the OWSC performance. Resonance of the sloshing modes is 
therefore beneficial for wave energy extraction. Despite the fact that the dynamics of the 
OWSC in the channel is affected by the bounding effect of the lateral walls, the general 
behaviour of the device in the open ocean can also be inferred. There the dominant 
physical mechanism is still the transmission of the incident waves beyond the flap, but 
no sloshing modes are present. Hence the OWSC is expected to perform in the open 
ocean similarly as it does in the 3D channel, exception made for the narrow bandwidths 
close to the resonant frequencies of the sloshing modes. However, to simulate properly the 
behaviour of a single device in the open ocean a different mathematical model is needed, 
in which the waves are allowed to propagate in all directions. In summary, care should 
be taken when using experimental results obtained in wave tank testing as a benchmark 
for predicting real on-site behaviour. In the following we shall investigate the influence 
of the system main parameter, namely the width ratio w = w' /b', on the performance of 
the OWSC. 

3.4. Parametric analysis and further research directions 

The response of the capture factor Cf (|2.41|) to variations of the flap width is shown in 
figure |U Here the plots of Cf against the wave period T 1 are shown for different widths 
of the OWSC, i.e. id'=6m (corresponding to w = w' jb' = 0.07), w' — 12m (w = 0.13) 
and w' = 18m (w = 0.2), for fixed width of the channel, b' — 91.6m. Note that for 
all widths the capture factor is l argest in short periods, thus matching the experimental 
observations of lWhittaker et al. ( 2007 ) . Indeed at larger periods (not shown here) , tuning 



of the incoming wave frequency with the natural pitching frequency of the flaps (|2.39b ) 
can occur, thus yielding additional peaks in the capture factor curves. U nfortunately, 



this e vent is unlikely to occur in normal operating conditions of the device (jHenrv et 



201Cft and won't be investigated here. Further examination of figure |S] also reveals that 



the wider the flap the larger the capture factor, for a g iven width of the c hannel. This 
funda mental result, already noticed experimentally by iHenrvl t|2008h and iHenrv et al 



(|2010h . reveals a selective behaviour of the sloshing modes with respect to the width of 



the flap, for which the most powerful resonance occurs with the largest flaps. Now, to 
what extent can we increase the flap width to exploit this mechanism, in a channel of 
fixed width? Of course, growth in the structural loads on the device is a heavy limit on the 
enlargement of the flap (jHenr v 2008). Hence we consider the alternative way of decreasing 



the width of the channel, while holding w' fixed. By decreasing b' so that w — w'/b' —> 1, 
the width of the channel becomes comparable to the width of the flap. Hence the 3D 
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Figure 8. Behaviour of the capture factor CV (|2.41|) versus the wave period T' for various flap 
widths. The solid line refers to w = 0.2 (w' = 18m in physical variables), the dashed line to 
w = 0.13 [w 1 = 12 m) and the dotted line to w = 0.07 (V = 6 m). 
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Figure 9. Behaviour of the excitation torque with the wave period in physical variables. The 
solid line shows the results for the 3D model with b' — 18 m and w' = 17.8 m, i.e. w — 0.99. The 
dashed line shows the results for the 2D model (|3 . 1 f) where b' — w' — 18 m, i.e. w = 1. 



effects become less noticeable, the geometry approaching that of a 2D device. Figure 
[9] shows a preliminary analysis of the torque F' with the incident wave period T' for 
w = 0.99. Here 3D effect s are largely inhibited b y the body, so that the device shows 
a "quasi-2D" behaviour (|Yeung fc Sphaierlll989h . As a consequence, when w = 0{\) 
the selective behaviour of the sloshing modes is expected to have a weaker impact on 
the performance of the OWSC, whose maximum efficiency will decrease closer to the 
theoretical 2D values. This suggests that an optimal layout maximising the resonant 
actions on the device must lie in between the two extremes w — > (open ocean) and 
w —> 1 (2D channel). A supplementary investigation of the response of the system to 



18 



E. Renzi and F. Dias 



variations of the channel width is then needed. This analysis is currently ongoing and 
will be disclosed in the near future. 

Finally, recall that resonance of sloshing modes is a linear mechanism. Nevertheless, 
resonant mechanism s can also excite the second-order transverse modes (cross waves, see 
IShemer et al. 1987 ). A preliminary numerical analysis with a weakly nonlinear Boussi- 
nesq model showed nonlinear amplification of the dynamic pressur e on the flap to occur 
near the cross- wave resonant frequencies (see ISchmitt et al\ 120 12h . A nonlinear mathe- 
matical model based on the governing equations of Appendix |A] and on the linear results 
obtained here is to be developed to investigate whether the nonlinear interactions can 
further increase the performance of the device. 



4. Conclusion 

A mathematical 3D model of an oscillating wave energy converter (i.e. the OWSC) 
in a channel is developed to investigate the influence of the channel sloshing modes on 
the performance of the system, in a layout common to m any experimental studies (see 
Follev et aZ]l2005t IWhittaker et aZ.ll2007t iHenrv et alJEnch . The linearised theory, based 



on the solution of a hypersingular problem in the fluid domain, is satisfa ctorily validated 
by comparis on with numerical and experimental models (see respectivelv lvan't Hofjl2009t 
lHenryll2008h . Whe n the flap is set to oscillations by incident waves, the incoming 2D wave 
field is modified by the diffraction occurring at the edges of the flap and by the waves ra- 
diating from the device. Then a complex 3D wave field develops in the channel, increasing 
the average efficiency of the OWSC with respect to that of a 2D pitching flap. In addi- 
tion, this 3D dynamics leads to the excitation of transverse sloshing waves which resonate 
each at its own cut-off frequency. Transverse modes trapped near the flap are responsible 
for the spiky behaviour of the hydrodynamic characteristics of the system and further 
increase the efficiency of the device near resonance. Finally, parametric analysis reveals a 
selective behaviour of the sloshing modes, whose resonant capacity depends on the width 
of the device. In response to the three issues we anticipated in the introductory section, 
we conclude that (i) resonance of the channel sloshing modes is a linear mechanism and 
occurs for every width of the flap; (ii) the sloshing modes show a selective behaviour 
with respect to the width of the flap, that concurs in increasing the performance of the 
device near resonance; (iii) care should be taken when employing the results obtained 
in an experimental wave tank to predict the behaviour of the device in the open ocean, 
especially near resonant frequencies. 

This work was funded by Science Foundation Ireland (SFI) under the research project 
"High- end compu tational modelling for wave energy s ystems" . Num erical simulations of 
Ivan't Hofj (120091 ) and experimental data provided by iHenrvl ([20 08) in agreement with 



Aquamarine Power Ltd have been very useful for the validation of the model. Fruitful 
discussions with Prof. P. Sammarco, Dr K. Doherty and Dr S. Bourdicr arc kindly ac- 
knowledged. The Authors wish to thank the wave research group at Queen's University 
Belfast for having let them partake in the June 2011 experimental campaign. 



Appendix A. Nonlinear governing equations 

Consider the geometry of figure [1] representing the OWSC in the channel as described 
in SJ21 During the motion, the position of each side of the converter (made up by the 
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bottom foundation and the flap) is at x' = X , where 

X' ± (z't') = ( ±Q '' - h '<z'<-h' + d 

v ' ; \ ~{z' + h'- d) tan 0'± a'/ cos 6', -ti + d < z' < Z'± ' v ; 

In the latter, the + (— ) superscript denotes the right (left) side of the OWSC and Z' + 
(Z'~) is the free-surface elevation along the wetted line of the flap on the right (left) 
lateral surface (see figure [TJi). Clearly, in (|A II) the upper equivalence holds for all points 
on the fixed bottom foundation, while the lower one is valid on th e lateral surface s of 
the moving flap. As in many gravity-wave problems (see for example Mei et aZll2Q05l ) the 
fluid is deemed to be inviscid and incompressible and the flow irrotational. Hence there 
exists a potential $'(x', y' , z', t 1 ) for the velocity field v' — V'$' such that 

V' 2 $' = 0, (x', y', z') G ft', (A 2) 

where f2' is the fluid domain. To set up a boundary- value problem governing the behaviour 
of the fluid in il', the Laplace equation (|A 2[) must be supplied with appropriate boundary 
conditions. On the free-surface, the kinematic-dynamic boundary condition reads 

&t, t , + g& M , + | V'$'| 2 t , + iV'$' • V | V'$'| 2 = 0, z' = C, (A 3) 

where g is the acceleration due to gravity. In (|A3[) the free-surface elevation £' is obtained 
via the Bernoulli equation 

-£ =92 < + *>, +i |V<*f, (A4) 

where p is the water density and p' \x' ,y' , z' ,t') the fluid pressure. By evaluating (|A 4[) 
at z' = £' and considering only the excess pressure over the atmospheric value, the 
free-surface elevation is 

C' = --4>: t * |V'<i>'| 2 , Z ' = C '. (A 5) 

Boundary conditions must also be applied at the solid frontiers delimiting the domain. 
We require absence of normal flux through the channel walls 

&y,=Q, y' = ±b'/2, (A 6) 

through the moving sides of the OWSC parallel to the x'-axis 

$; v =0, X'-(z',t') <x' < X'+(z',t'),y' = ±w'/2, (A 7) 

and through the impermeable bottom of the channel 

& z , =0, z' = -h'. (A 8) 

Finally, a kinematic boundary condition allowing only tangential motion alon g the lateral 
surfac es of the OWSC x' — X J± is to be applied. Following the reasoning of Sammarcol 
( 1996( ) for the mobile gates of the Venice storm barriers, we require d/dt' [x' — X /=t ] = 
[x - X'^lr + V'$' • V[x' - X f± ] = 0, which yields 

= J -4— [- ( z ' + h'- d) ± a' sinfi'] - tan0' \ 
\ cos 2 6' J 

xH(z' + h'-d), x' = X'±, ~ <y'<y, (A9) 

where (|A ip has been employed. In (|A 91) usage of the Heaviside step function H assures 
absence of flux through the bottom foundation. 
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The equation of motion of the flap 

I'e[ t , t ,=T' g {t')+T^t')+r{t') (A 10) 

expresses the dynamic equilibrium of torque about the hinge. In (|A 10[) I 1 is the second 
moment of inertia of the flap, assumed to be given, while T' g and T' are respectively 
the torque due to gravity and to the fluid pressure. Finally T' is the torque exerted on 
the flap by t he generator, wh ich in general is partly inertial, partly elastic and partly 
damping fsee lMei et al\ 120051 ) . Torque is positive when it makes the flap rotate counter- 
clockwise, according to the convention adopted for & (see figure Q}i). Since the fluid 
is considered inviscid and the hinge frictionless, there is no torque induced by viscous 
tangential stresses. The torque due to gravity is 

T' g {t') = S'gsin9'(t'), (All) 

where S' = M'd! is the first moment of inertia of the plate. In the latter expression M' is 
the mass of the flap and df is th e distance between the centre of mass and the hinge (see 
Mei et oill993 ls ammarc 1996). Both M' and d 1 are assumed to be given. The net torque 



T' v exerted by the fluid pressure p' on the flap is obtained by integrating the product 
between the unit pressure force p'dy'dz'/ cos 9' and the arm (z! + h! — c'^fa' sin 8') / cos 9' 
respectively on the right and left surfaces of the flap: 

rpi f Z ' + f w ' /2 / , ,a z' + h'-d-a'smQ' A fAj 

Z f"V(A" (A12) 

' —h'+c' J —iv'/l cos 

Summarizing, the complete set of nonlinear differential equations governing the coupled 
motion of the water and the flap includes: the Laplace equation (|A 2[) . the kinematic- 
dynamic boundary condition (|A 3j) at the free surface (|A 5| . the no- flux conditions re- 
spectively on the channel walls (|A 6j) . on the sides of the OWSC parallel to x' (|A 7^ and 
on the bottom of the channel (|A 8|) . the kinematic condition on the lateral surfaces of 
the OWSC (|A9j) and finally the equation of motion of the flap (|A 10[) . 



Appendix B. Solution of the plane radiation and scattering problems 

Solution to the plane radiation and scattering problems of £12.21 is sought here with an 
integral-equation technique based on an appropriate decomposition of the channel Green 
function. 

B.l. The channel Green function 

Define the Green function G n (x,y;!;,r]) singular at (£,±77) € S as the outgoing solution 
of the Helmholtz equation 

(V 2 + K 2) Gn = , (x,y) e£\(£,±r?), (B 1) 

where S is the 2D fluid domain of figure Wlb), with boundary conditions 

G n , y =0, y = ±l/2, (B2) 

G„~-^lnp ± , p ± ^0, (B3) 
Ztt 

where = \J {x — £) 2 + (y T v) 2 an d in are the eigenvalues ()2.23j) of the dispersion 
relationship (|2.24[) . n = 0, 1,2, . . . The solution of the boundary-value problem (|B lj) ~ 
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(IB 3j) is given in its classical form by iLintonl ( 19981 ) . Here we shall use the decomposition 



G n (x, y; £,??) = — 7 [g£ 0) [x, y; £ r,) + G« (x, y; £, rf) 



(B4) 



where 



Gi°\x,y;tv) = («„ y/(x - £) 2 + (v ~ v) 2 ) + ^ («n V 7 ^ - £) 2 + (y + ^) 2 ) 



(B5) 



and 



G«(x, V ;f^)= [ H o 1] (*nV(x-0 2 + (y-V-™) 2 ) 

X) 

I 



m— — oo 
m/0 



(B6) 



with Hji' the Hankel function of first kind and order n € N. Expressions (IB 4I) - (|B 6p 
represent a sum of sources equally spaced along the y-axis, due to the mirroring effect of 
the side walls (as anticipated in q2.1[) . Note that Gn (|B 51) is singular at (x, y) — (£, ±n) G 
E, while G^ 1 (|B 6|) has no poles in E. The simple Green function decomposition (|B 4j) into 
a singular part and a converging series will enable us to find appealing fast-convergent 
semi-analytical solutions to the radiation and diffraction problems. 



B.2. Solution 

Consider the plane potentials ' D ' solving the system (|2.25p ~ (|2.27p . outgoing at large 
|at|. Application of the Green theorem (see for example Meil 119971 ) to tfn 1 ' ^ and G n in 
the domain E yields after some algebra 



v/2 



w/2 



G 



(0) 



G 



(i) 



dr?, 



where 



A^> = A^>fo) = ^(-O,,) - ^ D \+0,v) 



(B7) 



(BJ 



denotes the jumps in radiation and diffraction potentials from the left to the right side of 
the flap in the x direction, still unknown. In (|B 7[) the subscript £ = indicates the point 
at which the term in brackets is to be calculated. Application of the boundary conditions 




G 



(0) 



G 



(i) 



dr) 



gives 






f ^/« 1 


x=±0 


_ A <2„ J 



-w/2 < y < w/2, 



(B9) 

which are integro-differential equations for Aip n R ' D \ The structure of (|B 9[) would sug- 
gest to bring the outer derivative under the integral sign. However, due to the singular 
behaviour of G^ (IB 51). this would l ead to a divergent integrand near the poles r\ = ±y. 



Nevertheless, Martin fe Rizzol ( 19891 ) demonstrated that in integro-differential equations 



like (|B 9| the inversion between the outer derivative and the int egral is possible provide d 
the latter is interpreted as a Hadamard finite-part integral ^ (jLinton fc Mclver 2001 ). 
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By virtue of Martin fe Rizzo ( 1989t l's theorem, (IB 91) can be rewritten as 




(BIO) 



The latter are hypersingular integral equations for the jumps in potential Atp^' D ^ across 
the flap. The singularity lies in the first term of the left-hand side and is due indeed to 
the divergent component G n °^ of the Green function G n (|B 4[) . To reveal the nature of 



the singularity, substitute the series (lB5l) and ([B6j) respectively for G ( „ ' and G n ' into 
(|B 10p . Then perform the double differen tiation with respect to £ an d x and use the 



property /" q f(x) dx = f° a f(-x) dx (see iGradshtevn fe Rvzhikll2007l §3.022) to get 



P n (u) \ (±K n w\v - u\) 



+ 00 

E 

m— — oo 



Pn(u) 
Qn(u) 



Qn(u) 

Hf' (±K n w\v m - 



du 



u\) 



du 



\vo - 
4i 



A d n 



for the boundary condition on the wet contour of the flap. In the latter, 

= 2(y-m)/w, r),y G (-w/2,w/2), 



while 




(Bll) 



(B12) 



A^(r/) 



(B13) 



Atp^wu) 1 

denote the jumps in potential in the new variables. Note that P n (u) = P n (—u) and 
Q n {u) = Qn{—u) because of symmetry. Again, as a result of the decomposition (|B 4|) . 
the singular behaviour of (|B 111) is restricted to the finite-part integral of the left-hand 
side. We must focus on this term to circumvent the singularity. First, expand the Hankel 
function according to the series representation (see IGradshtevn fc Rvzhikl [20071 

§8.444) 

4 1 

R n (\k„w\v - u\) , 



H^\±K n w\v Q - u\) = 



ITT K n W\V() — U\ 



(B14) 



where 



R n (a) = Ji(a) 



1 



7T V 2 



+ 7 



2 + ^ 

i=2 



J/'Ki-l)! 



E 

9=1 



9 



(H 15) 

is the remainder, Ji(x) is the Bessel function of first kind and first order and 7 = 
0.577215 . . . the Euler constant. Then, by making use of (IB 141) . rewrite the hypersingular 
integral equations (|B lip as 



Pn{u) 
Qn{u) 



(vq — u) 2 du + 



IC n [Pn(u)} 
IC n [Qn(u)} 



Vfn 
A d n 



(B16) 
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In the latter expression 



ICn[f(u)} = 



/(«) 



R n (±k n w\v - u\) 



\vo 



E 



m^O 



du 

(B17) 

is an integral function. Note that as u — > vq G (—1,1), R n (^K n w\vo — u\) ~ \vq — 
w| In |wo — u| from (|B 15|) . hence IC n has a convergent kernel. The singularity in (|B 16)) is 
finally isolated in the Hadamard integral, whose kernel has a 2nd-order pole at vq. This 
fundamental result now allows us to solve the hypersingular integral equation (|B 16|) 
in terms of the jumps in radiation and diffraction potentials, P n and Q n respectively. 
Indeed the Hadamard integral in (|B 16[) admits eige nfunctions that are pro portional to 
the second-kind Chebyshev polynomials U p , p G N ( Linton fc Mclverl 2001 ). The latter 
constitute a complete set of eigenfunctions in the domain u G [—1,1]. Then the analytical 
form of (IB 16|) itself suggests to seek for solutions of the type 



Pn{u) 
Qn(u) 



V 

A 



1-u 



+ 00 

- 1/2 E 

p=0 



Qpn 
fipn 



U p {u), 



(B18) 



where the a pn and /3 pn are unknown complex constants. Note that the position (|B 181) 
also reflects the asymptotic behaviour of the jump in potential near the tips of the flap 
u = ±1, where P n , Q n ~ \/ (1 — u 2 ) — > an d the velocity has a square-root singularity 



(see iNewmannll 1 9 7 It iLinton fc Mclver 200 l|). Now the Chebyshev polynomials U p satisfy 
the integral relationship (see lParsons fc Martini 19921 ) 



u 2 ) 1/2 Up(u) 



du = -7r(p + l)Up(v ), w 6 (-1,1), 



(B19) 



where vo = 2y/w according to (|B 12[) and y E (— w/2, w/2). Substitution of the series 
expansions (|B 18)) in the relevant hypersingular integral equations (|B 16)) and application 
of the property (|B 19)) yield finally 



oo 

E 

p=0 



Cpn(v ) = -irw 



with 



C pn (v ) 



-ir(p+l)Up(v Q ) + 

R n (iKnW\v ~ u\) 

| wo - u\ 



In 



(i 



e (-1,1), 



(B 20) 



u 2 ) 1/2 U p (u) 



E 

m— — oo 



H± (^K n w\vo — 2m/w — u\) 
| Wo — 2m/w — u\ 



du. (B21) 



The latter can be evaluated by solving numerically the convergent integral in (|B 21)) . 
For given n, expression ()B 201) defines two different linear equations, both valid for any 
wo € (—1,1). The upper equation is for the coefficients a pn = a pn (uj) of the jump in 
the plane radiation potential ip^ across the plate, while the lower equation is for the 
coefficients f3 pn = /3 p „(u;) of the jump in the plane diffraction potential (p„. Theoretically 
there exists an infinity of points wo G (—1, 1) where (|B 20[) can be evaluated, thus yielding 
a system of an infinite number of equations (one for each wo) for an infinite number of 
unknowns. For (|B 20)) to be solved numerically, the unknowns must be truncated to 
a finite integer number P < oo and a finite number of evaluation points wq = wqj, 



24 
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j = 0,1,..., P must be also chosen in the domain (—1,1). Parsons fc Martini (1992) 



3 

showed that the fastest numerical convergence for such a syste m is achieved when the 



voi a re the zeros of the Chebyshev polynomial of first kind (see iGradshtevn fc Rvzhik 



2007), i.e 



(2i 4- 1W 

v V =cos 2P + 2 ' 3 = 0,1,..., P. (B22) 
Hence (|B 20j) reduces to two (P + 1) x (P + 1) truncated algebraic systems 

E{ ^ }cWt*w) = -™{ £ }, J=0,1,...,P, (B23) 

for each n = 0, 1, 2, . . . Once (IB 23[) has been solved in terms of the a pn and /3 pn for a 
given modal order n, the jump in the ra-th modal potentials across the plate Aipif'' ^ (rj) 
can be determined with (|B 18[) together with the change of variable dictated by (|B 12|) 
and (|B 13|) . Now substituting Aip^' ^^) into Green's theorem (|B 7[) . differentiating the 
Green function (|B 4j) inside the integral and summing up all the eigenmodes yield the 
expressions of the plane potentials. The n-th modal plane radiation potential is given by 

<Pn( x iV) = ^~K n X ^a (2p )„ ^ / (l-U 2 ) U 2 p(u) 

p—0 m— — oo — ^ 

ij} 1 ^ ( K n ^J X 2 + (y — — m) 2 ) 

x \ „ . ■ — = L du, (B24) 

yx + \y — ±wu — my 

where M — [P/2\ . In (|B 24[) only the even terms survive because the odd terms give 
no contribution when integrated and summed up over symmetric domains. Finally, the 
plane diffraction potential is given by 

tfio(x,y)= kx ^2 /3 {2p )o / (1 - u 2 ) 1 2 U 2p (u) 

p—0 m— — oo — 1 

h[ (kyjx 2 + (y — \wu — m) 2 \ 

x V , 2 = ■ = ' du. (B25) 

+ (y — jWM — m) z 

Here only the 0-th modal component tp® is non-zero; in other words, the diffraction 
problem only admits the fundamental mode n = 0. This happens since d n = for n > 
(see 12.301) . Hence the /3 pn , which are homogeneous solutions of (|B 23[) for n > 0, must 
be zero to ensure uniqueness of the non-homogeneous solutions a pn for n > 0. This is 
a solvability condition of the coupled radiation-diffraction problems, which substituted 
into (IB 18[) and then into (|B 7p yields ultimately ip„ = 0, n > 0. 



Appendix C. Further relations between radiation and scattering 

In this section we provide some relations between the radiation and the scattering 
problems discussed in Sj^l First consider the excitation torque F given by (I2.36|) . The 
latter depends on /3qo, solution of the linear system (|B 23|) for the scattering problem. 
Here we show that that F can also be written in terms of the solutions «oo to the linear 
system (|B 23|) for the radiation problem. Consider the two systems of equations (|B 23j) for 
a p o and /3 p o respectively. Since the coefficients C p o are the same in either system, then 



Resonant behaviour of an OWEC in a channel 



25 



one must have /3 p o = ct p odo //o to ensure uniqueness of the solution. Hence substitution 
of the latter expression into (|2.36|) and the use of (|2.30p for do yield after some algebra 

i7T«; Vk{2kh + smh2kh) 1/2 
F --— Ao c^shkh ao ° (C1) 

for the complex excitation torque. In (jC II) aoo is given by the solution of the system (|B 23|) 
for the radiation problem, which then suffices to describe the dynamic actions on the 
flap. In addition, a direct relationship between the diffraction torque F and the radiation 
damping v (|2.35|) can also be found. First, numerical evaluation of the coefficients ao„ 
reveals that Imjaon} = for all n > 0. As a consequence, expression (|2.35l) for the 
radiation damping simplifies as 

^=™/oIm{a o}. (C2) 

Then by taking the real part of F (jC 1[) and making use of (jC 21) to express Im{aoo} in 
terms of is, the sought relationship between the radiation damping and the excitation 
torque is found 

cosh kc — cosh kh 



1 



k Cg 



c + 



k sinh kh 



hm]i/,7/ Rc<{ \ . (<?:-!) 



C g being the non-dimensional group velocity corresponding to (|2.42j) . Expression (|C 3j) is 
the counterpart of the 2D Haskind-Hanaoka relation that associates the excitation torque 
acting on the hxed flap in incoming waves with the r adiation potential of the pitching 



motion of the body (see for example I Mei et all 120051 ) . According to (jC 3[) a close link 



exists between the radiation and the scattering problems, for which they both exhibit 
the same resonant behaviour, as pointed out in 8 33.21 Finally, (jC 3| has also been employed 
as a benchmark to double-check the solution of the whole system. 
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